Temperature dependence of the thermal boundary resistivity of glass-embedded metal 

nanoparticles 
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The temperature dependence of the thermal boundary resistivity is investigated in glass-embedded 
Ag particles of radius 4.5 nm, in the temperature range from 300 to 70 K, using all-optical time- 
resolved nanocalorimetry. The present results provide a benchmark for theories aiming at explaining 
the thermal boundary resistivity at the interface between metal nanoparticles and their environment, 
a topic of great relevance when tailoring thermal energy delivery from nanoparticles as for applica- 
tions in nanomedicine and thermal management at the nanoscale. 
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With the ever decreasing size of nanodevices, investiga- 
tion and modeling of heat exchange at the nanoscale has 
become of central technological interest. Under a fun- 
damental standpoint, metal nanoparticles (NP) embed- 
ded in a host matrix constitute a model system, as they 
can be selectively heated-up and their cooling monitored 
using time-resolved spectroscopy^—. Furthermore, the 
thermal dynamics occurring between an optically excited 
metal nanoparticle and the surrounding environment is 
of direct relevance for a variety of applications rang- 
ing from photothermal cancer therapy 4 ^ and selective 
drug delivery^ to thermoacoustic imaging and electro- 
magnetic waveguiding in dielectric-embedded plasmonic 
devices J- The electromagnetic energy harvested by the 
NP is dissipated as thermal energy in the environment. 
The corresponding energy flux J p is ruled by the thermal 
boundary resistivity pf,d, i-e., Kapitza resistivity, and by 
the temperature mismatch AT between the two media: 
J P =AT / pbd- Investigating pm is therefore crucial to tai- 
lor thermal energy delivery from the NP to the matrix or 
matrix-embedded target and, more generally, to analyze 
heat transfer at the nanoscale. 

Whereas effort has been devoted to understand and 
model the Kapitza resistivity between two solids, both 
bulk and thin films^— , the scenario remains relatively 
unexplored when one of the two materials downscales to 
the nanometer range. Several confinement effects may 
modify p^d, for instance, as the dimension of the NP 
becomes comparable to the thermal diffusion length of 
the host material 1 - 1 -"— , or is reduced to a point where 
the continuum solid approximation to the elastic problem 
becomes questionable. A lack of extensive experimental 
evidence^ ' 14 i 15 spanning the space of parameters affect- 
ing p^, most notably the temperature^, has so far pre- 
vented a consistent account of the mechanisms ruling the 
Kapitza resistivity at the nanoscale. When confronted 
with the problem of measuring the heat transfer from a 
nanoscale object, the challenges stand in: (a) a probe 
speed requirement, dictated by the fact that the time for 
heat exchange between the sample and the thermal reser- 
voir decreases with the decreasing sample's mass; (b) a 
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FIG. 1: (Color online) Measured time-resolved relative trans- 
mission change for T cry o = 200 K. The pump and probe pulses 
have wavelengths respectively at 400 nm and 800 nm. In 
the cartoon the thermal fluxes J p and J m are represented to- 
gether with the temperature profile within the sample. Inset: 
measured OD of the sample outlining the Ag nanoparticles' 
LSPR. 



non-contact probe requirement, to avoid the addendum 
heat capacitance contribution from the probe itself. 

In this Letter we use time-resolved all-optical 
nanocalorimetry^ to overcome such challenges and in- 
vestigate the temperature dependence of the cooling dy- 
namics of glass-embedded Ag particles of radius R — 4.5 
nm. The Kapitza resistivity is shown to increase by a 
factor of two with decreasing temperature from 300 to 
70 K, a trend consistent with existing models. 

The Ag nanospheres are embedded in a 50% BaO - 
50% P2O5 glass matrix of thickness L — 50 /im. The 
metal volume fraction is 2 • 10~ 4 . The sample was syn- 
thesized using a fusion and heat treatment techniquej 17 i 18 
The samples' optical density (OD) shows enhanced ab- 
sorption in the blue portion of the spectrum due to the 
localized surface plasmon resonance (LSPR) of the Ag 
NP - see inset of Fig. [TJ 

The time-resolved measurements were performed using 
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a Ti: Sapphire cavity dumped oscillator - 800 nm wave- 
length, 120 fs pulse temporal width at full width half 
maximum. The Ag NPs are selectively excited by the 
frequency-doubled pulse at 400 nm wavelength close to 
the LSPR in order to maximize energy absorption in the 
particle. This leads to a fast heating of the electrons of 
the NPs that thermalize with the lattice on a few picosec- 
onds time-scale, the thermal energy being subsequently 
delivered to the matrix. Care was taken to minimize av- 
erage heating-is of the glass matrix by keeping the energy 
per pulse as low as possible while granting a detectable 
transmission variation. To this end the laser repetition 
rate was tuned to 540 kHz by means of a cavity dumper. 
Accounting for transmission losses along the optical path, 
the energy density per pump pulse on the sample surface 
was Io ~ 0.5 J/m 2 , while the energy density absorbed per 
particle per pulse was U v ~ 4-10 7 J/m 3 . The cooling dy- 
namics of the hot NPs to the glass matrix is then followed 
by measuring the relative change in transmission across 
the sample ATr/Tr of a time-delayed 800 nm probe pulse. 
Probing out of the LSPR grants proportionality between 
the experimental signal and the NPs temperature rise^ 
(at the expense of the signal amplitude) . 

A typical experimental trace is reported in Fig. [T] for a 
cryostat temperature T cryo — 200 K. After excitation by 
the pump pulse and internal electron-lattice thermaliza- 
tion, i.e., after 6 picoseconds (this step has been exten- 
sively investigated in these systemai^ and will not be 
discussed here), the signal decay reflects cooling of the 
hot NP to the matrix. This is governed by the thermal 
flux J p and J m from the NP to the matrix and from the 
matrix portion adjacent to the NP to the rest of the ma- 
trix, respectively. Considering these two processes, the 
energy balance is governed by: 



C p dtT p (t) = — 



PX ' Rpbd 

C m d t T m {r,t) = A m r- x ^[rT m (r,i)] 



[T p (t)-T m (R,t)] (1) 
(2) 



T p being the NP's temperature, assumed as constant 
throughout the particle volume, T m the matrix' tempera- 
ture, Cp and C m the particle's and matrix's specific heat 
per unit volume respectively, and A m the matrix' thermal 
conductivity. For the case of constant thermal parame- 
ters, the temperature increase for the NP, AT p (t), and 
for the matrix portion in contact with it, AT m (R, t), are 
analytically accessible working in Laplace space 2 ^ and 
read: 
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2k(Rg) 2 AT 



u 2 exp(-KU 2 t/R 2 



[u 2 (l + Rg) 



kRg] 2 + (it 3 - kRgu) 2 
(5) 
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FIG. 2: (Color online) Time evolution of temperature and 
specific heat for T cryo = 70 K. Panel (a): relative temper- 
ature variation (left axis) and absolute temperature (right 
axis) of the NP (full line) and of the adjacent matrix (dashed 
line). Inset: experimental transmission change normalized to 
the value at 6 ps (red curve) and its best fit (black curve). 
Panel (b): relative specific heat variation (left axis) and abso- 
lute specific heat (right axis) of the NP (full line) and of the 
adjacent matrix (dashed line). 
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AT is the initial temperature increase 2 ^ and 
k = A m /C m , k = 3C m /Cp, g = l/A m p bd . In our 
analysis the thermal resistivity pbd is set as a fit param- 
eter, together with A m , which is not precisely known for 
our glass material. 

As low temperatures are investigated (in particu- 
lar, around and below the NP's Debye temperature, 
Td ~ 215 K for Ag), C p and C„r& cannot be set to 
their T cryo value and regarded as constant over the 
particle and matrix's temperature excursion taking place 
during the experiment. The solution of Eq.s [1] and [21 
taking into account the temperature dependent specific 
heats, is then retrieved iteratively. Fitting is performed 



starting at a time-delay of 6 ps and setting T, 



m,0 



T, 



and T, 



cryo 



f ATq. The corresponding values 
G p {Tpfi) and C m (T mi o) are inserted in Eq.s [3] and |4j 
and the new temperatures T Pj i = T cryo + AT Pj i and 
T m ,i = T cryo + AT m .i are calculated and adopted in 
the subsequent time step. The procedure is iterated 
to reach the maximum experimental time-delay of 320 
ps. Values for p^d and A m are obtained maximizing the 
likelihood between the theoretical T p (t; pbd, Am)/ATo 
and experimental —ATr/Tr traces. 

The resulting dynamics of temperatures and specific 
heats are exemplified in Fig. [2] (a) and (b) respectively 
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FIG. 3: (Color online) Kapitza resistivity pbd vs Tcryo (black 
circles). The horizontal arrows indicate the temperatures 
spanned by the NP during the thermalization process. Plot 
of the function ACp 1 , A being a multiplication constant with 
dimensions ms _1 (red curve). Inset: normalized transmission 
change (red curve) and its best fit (black curve) for the case 
T cryo =200 K. 



for the lowest studied temperature. The internal 
thermalization of the NP is achieved at T p — 93 K. As 
time evolves the NP cools down, increasing T m (R), the 
maximum value of T m (R) being attained at t ~ 15 ps. 
For longer time-delays both T p and T m (R) decay toward 
the asymptotic value T cryo . The same trend applies to 
the specific heats C p and C m (R), showing maximum 
relative changes during the experiment in the 15-20% 
range. A similar behavior is obtained for measurements 
performed at higher temperatures, although with a 
smaller excursion amplitude (as indicated by the arrows 
in Fig. [3]). These variations stress the importance of 
taking into account the temperature dependence of the 
specific heat when measuring at cryogenic temperatures 
and make difficult the extraction of ptd at values of 
Tcryo < 70 K. The Kapitza resistivity Pbd(T cryo ) 



increases by a factor of two, spanning values from 
3.2 to 6.5 m 2 K/GW, as the cryostat's temperature 
decreases from 300 K to 77 K, see Fig. [3] The extracted 
value corresponds to a mean value of p^d over the 
NP's temperature excursion during the experiment 
(indicated by arrows in Fig. [3]). Values for A m were 
found in the range 0.2-0.7 W/mK, comparable to the 
ones reported for thermal conductivity of glasses with 
similar compositions^ 

Starting from the general expression for the Kapitza 
resistivity-^, and assuming both a frequency independent 
and/or frequency-averaged phonon transmission coeffi- 
cient t, and group velocity v g , one finds pm ~ {tv g C p )~ l . 
This trend is experimentally retrieved in our data 
where ptd is found to roughly follow the temperature 
dependence of C" 1 , see Fig. [31 

In conclusion we measured via time-resolved all- 
optical nanocalorimctry the Kapitza resistivity of a 
4.5 nm radius glass-embedded Ag nanoparticle in the 
temperature range from 300 K to 70 K. p^d increases 
monotonically from an ambient temperature value of 
3.2 m 2 K/GW to 6.5 m 2 K/GW. The present findings 
constitute a benchmark for theories aiming at explaining 
the Kapitza resistivity in nanosystems, a fundamental 
issue for applications in nanomedicine and thermal 
management at the nanoscale. 
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